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Abstract. We verify a conjecture on the structure of higher-rank 
numerical ranges for a wide class of unitary and normal matrices. 
Using analytic and geometric techniques, we show precisely how 
the higher-rank numerical ranges for a generic unitary matrix are 
^f) ] given by complex polygons determined by the spectral structure 

of the matrix. We discuss applications of the results to quantum 
(N ■ error correction, specifically to the problem of identification and 

^ ' construction of codes for binary unitary noise models. 

(N 
oo 

^ ! 1. Introduction 

p . 

The study of higher-rank numerical ranges of matrices was initiated 
Oh! in [6], with a basic problem in quantum error correction [5] giving 

^ I the primary motivation. Higher-rank numerical ranges generalize the 

■ classical numerical range of a matrix, and arise as a special case of the 
matricial range for a matrix [[8l I15| . In ||6], three of us conjectured that 
the higher-rank numerical ranges of normal matrices depend, in a very 
precise way, on the spectral structure of the matrix. The conjecture 

^ ' reduces in the rank-1 case to a well-known property of the classical nu- 

■ merical range for a normal matrix, and it has opened the door to some 
interesting new mathematical problems. Its verification (or refutation) 
would also yield information for quantum error correction. 

In this paper we verify the higher-rank numerical range conjecture 
for a wide variety of unitary and normal matrices. This is accomplished 
by introducing a number of new geometric techniques into the analysis. 
We show in the case of a generic N x N unitary with non-degenerate 
spectrum and positive integer k > 1 with > 3k, that the kth numeri- 
cal range is given by a certain polygon in the complex plane determined 
by the eigenvalues of the unitary, and thus verify the conjecture. In 
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other cases, such as = 5m and k = 2m, we also verify the conjec- 
ture. But the analysis in these cases is more delicate, and our proof is 
non-constructive in nature. Figure 1 below provides a chart indicating 
the cases we verify, together with the various constraints. 

Our results may be applied to construct error correcting codes for a 
special class of quantum channels. A "binary unitary channel" [5] S is 
a noise model described by two unitary errors that can occur during the 
time evolution of a given quantum system. The construction of codes 
for such a channel relies on the structure of the higher-rank numerical 
ranges for a single unitary U. Suppose U acts on A^-dimensional Hilbert 
space. Then an "[A^, fc]-code" for £^ is a fc-dimensional subspace code 
that is correctable for S. The results described above for N > 3k yield 
a simple algorithm to determine the existence of codes, and an explicit 
construction of codes when they exist. 

The paper is organized as follows. Section 2 contains a discussion of 
some basics of quantum error correction that give important motivation 
for consideration of higher-rank numerical ranges. In Section 3 we 
discuss the conjecture, and show in particular how the general normal 
case relies on the unitary case. Section 4 deals with the structure 
of the aforementioned polygon and a derivation of conditions under 
which it is nonempty. Section 5 includes the proof for N > 3k and 
the construction of codes for binary unitary channels. In Section 6 we 
derive a number of other cases non-constructively based on the case 
N = 5 and k = 2. We finish with a brief discussion on possible further 
avenues of research and limitations in Section 7. 

2. Error Correction in Quantum Computing and Binary 

Unitary Channels 

For a more complete discussion on the material of this section see [5] 
and the references therein. We start with a quantum system in contact 
with an external environment having finitely many degrees of freedom 
represented on a Hilbert space H such that dimTi = A^ < oo. Consider 
a unitary time evolution of the combined system and environment in- 
duced by a given Hamiltonian associated with some quantum computa- 
tion implemented on Ti. The action of the evolution map on the system 
is obtained by tracing out the environment, and the resulting map is 
called a quantum channel or operation. Such a map is described by a 
completely positive, trace preserving map £ : C(H) — > C(TC), and can 
always be represented in the operator-sum form as £{p) = Ylia^apE\ 
for a set of operators {Eg} C £(7Y) satisfying '^a.-^a.Ea = I- As a 
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convenience we shall write £ = {Ea} when the operators Ea determine 
S in this way. The Ea are interpreted as the noise or errors induced by 
£. 

In the standard approach to quantum error correction, a code on Ti 
is given by a subspace C C 7i of dimension at least two. Denote the 
projection of Ti. onto C by Pc- A code C is (ideally) correctable for an 
operation £ if there is a quantum operation 71 : C{T-C) — > >C(7i) that 
acts as a left inverse of S on C; 

(1) {nos){p) = p wpePcC{n)Pc. 

Given a representation for S = {Ea}, a code C is correctable for S if 
and only if there are complex numbers {Xab) such that 

(2) PcElEhPc = XabPc ya,b. 

Thus, the problem of finding correctable codes for £ is equivalent to 
simultaneously solving the family of equations in Eqs. ([2]) for the scalars 
Xab and projections Pc, for all pairs a, b. 

Of central importance in quantum computing, communication and 
cryptography is the class of randomized unitary channels |JLJ, i3j. Such 
a channel has a representation of the form S = {^/fhiUa} where each 
Ua is a unitary operator and the {pa} form a classical probability dis- 
tribution; Pa > 0, J2aPa = 1- Heuce, 

(3) £ip) = J2paUapui ypecin). 

a 

The associated quantum operation is given by the scenario in which 
the error Ua occurs with probability pa- By Eqs. (|2]), finding ideal 
correctable codes for 8 = {^/p^Ua} is equivalent to solving the (un- 
normalized) equations 

(4) PUlUbP = XabP ya,b, 

for Xab and P. Note that each operator f/]f/fe is unitary. 

The following observation illustrates the importance of randomized 
unitary channels in error correction. 

Proposition 2.1. Let S be a quantum operation. Then C is a cor- 
rectable code for S if and only if there is a randomized unitary channel 
T such that C is correctable for T and 

(5) E{p)=J={p) \lpePcL{:K)Pc. 

Proof. In fact, T = {^/p^Ua} can be chosen so that the partial isome- 
trics UaPc have mutually orthogonal ranges for distinct a. This follows 
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directly from the usual construction of a correction operation TZ for 
£ on C |10] . The matrix {\ab) from Eqs.([5]) is a density matrix, and 
the unitary which diagonalizes it can be used to find a set of error 
operators {Fa} that implement £ o P^, where Pc(-) = Pc{.-)Pci such 
that PcF^FjjPc = SabdaaPc- The polar decomposition yields a partial 
isometry Ua, which can be extended to a unitary on the entire space, 

such that FaPc = UaV PcFlFaPc = \fdZUaPe- The form a 

probability distribution by the trace preservation of ■ 

Of course, the randomized unitary channel given by the restriction 
of £ to the code C can only be obtained when the code itself is known. 
Thus, this result in itself is of little practical use for the general prob- 
lem of finding error correcting codes. Nevertheless, it indicates that 
the general problem is equivalent to finding codes for the class of ran- 
domized unitary channels. 

Observe that a code is correctable for a binary unitary noise model 
£ = {^/pV, a/1 — pW}, where V and W are unitary, if and only if it is 
correctable for £ = {^/pl, \^1~^V^W}. 

Definition 2.2. A binary unitary channel on a Hilbert space Ti is 
a channel of the form 

(6) £ = {^I,./T^U} 

for some unitary U G C(H) and fixed probability < p < 1. Thus, the 
action of £ is given by 

(7) £ip)=pp+{l-p)UpU^ \/peC{n). 

Remark 2.3. Binary unitary channels form a rather restrictive class 
of physical noise maps, but they provide a useful set of "toy" examples 
for testing the "compression" approach [5J to build quantum error cor- 
recting codes enabled by consideration of higher-rank numerical ranges. 
Observe that from Eqs. the problem of finding ideal correctable 
codes for a given binary unitary channel is equivalent to solving four 
equations, but that the entire problem reduces to solving the single (un- 
normalized) equation for A and P given by: 

(8) PUP = XP. 

An immediate consequence of what follows is an algorithm to construct 
codes for a wide class of binary unitary channels. This is encapsulated 
in the discussion of Section 5. 
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3. Higher-Rank Numerical Range Conjecture 

Given a fixed positive integer k > 1 and T G C{T-C), the kth numerical 
range of T is the set of complex numbers 

A;.(T) = {A G C : PTP = XP for some rank— A; projection P}. 

The classical numerical range W{T) = Ai(T) is obtained when k = 1. 

Definition 3.1. Let T G C{T-C) and let k > 1 be a fixed positive inte- 
ger. Then we define Qk{T) to be the intersection of the convex hulls 
conv(r), where T runs through all {N — k + l)-point subsets (counting 
multiplicities) of the set of eigenvalues spec(T) for T. That is, 

^k{T) = fl conv(r). 

rCspec{T); |r|=Ar-fc+l 

Thus, f2fc(T) is a convex subset of the complex plane that can be 
computed directly from the spectrum of T. Below we will show how 
this set can typically be computed as an intersection of much fewer 
than sets. 

It is easy to see that flk(T) contains Afc(T) for normal T. We include 
a short proof for completeness and notational purposes. 

Proposition 3.2. Let T G C{T-C) be a normal operator and fix a positive 
integer k > 1. Then Ajt(T) C r2fc(T). 

Proof. Let {iV'i); ■ ■ ■ , \4'n)} be a complete set of orthonormal eigen- 
vectors for T with eigenvalues T\ipj) = Xj\ipj)- Let A G A^iT) and 
let P = J2i=i be a lank-k projection such that PTP = XP. 

Then {T(j)i\(j)j) = 6ijX for all i,j. Let A be a subset of {1, ... , A^} with 
cardinality \A\ = k — 1. Choose a unit vector |0) in the /c- dimensional 
subspace PH = span{|0i), . . . , \(pk)} that is perpendicular to all the 
lipj) for which j G A; and so, |0) = J2jfA^j\i'j) "^i^^ J2j = 1- 
Then we have 

(9) A=(T0|0) = 5^|z,fA, G conv{A, : J ^ A}, 

and it follows that A belongs to Qk(T). ■ 

It is well-known and easy to verify that the numerical range of a 
normal operator T coincides with the convex hull of its eigenvalues 
(that is, Ai(T) = Qi(T)). In [6J the following conjecture was asserted 
as a generalization of this fact. 
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Conjecture A. Let Ti be an A^-dimensional Hilbert space and let 
/c > 1 be a positive integer. Then for every normal operator T G C{T-C), 

(10) Afc(T) = Q,{T). 

The Hermitian case [B] and the normal < 4 case |5| of the conjec- 
ture have been verified previously. We show that the general normal 
case of Conjecture A can be reduced to the unitary case. 

Proposition 3.3. Conjecture A holds if and only if the conjecture holds 
for all unitary matrices. 

Proof. First note that for a fixed k, a standard translation argument 
shows the statement Afc(T) = Qk{T) for all normal T G C{C^) is 
equivalent to the statement G Afc(T) if and only if G ^k{T) for all 
normal T G C{C^). We focus on the latter formulation. 

Every normal operator T decomposes as T = Ti©Om, where Ti is nor- 
mal and invertible. The case m > k is easily handled, so assume m < k. 
One can check that G Afc(T) if and only if G Ak-m{Ti), and G 
flk(T) if and only if G Qk-m(Ti). Let {Ai, . . . , \N-m} be the (non- 
zero) eigenvalues for Ti, and let U be the unitary on C^""^ obtained 
from the polar decomposition of Ti with eigenvalues { • • • , \\^~"'\ }- 
By assumption we have Ak-m{U) = Qk-m{U), and hence zero belongs 
to both sets or neither set. Thus, we complete the proof by show- 
ing that: (i) G Ak-m(Ti) if and only if G Ak-m{U), and (ii) 
G flk-m{U) if and only if G Qk-m(Ti). 

Note that by invertibility we have Ti = UR = RU = \fR U VR 
where R = VT^T > is invertible with eigenvalues {|Ai|, . . . , |Ajfc_rrt|}. 
Thus, (i) follows from the more general principle that if T = X'^SX 
where X is invertible, then G Afc(T) if and only if G Ak{S). Indeed, 
if P is a rank-fc projection such that FTP = 0, then the rank-fc range 
projection Q of XP satisfies QSQ = 0. 

For (ii), note that by definition ^ Qk-miTi) precisely when does 
not belong to the convex hull ofA^ — fc + m+ lof the eigenvalues 
{Ai, . . . , Xn-tu}- This is equivalent to the existence of a line passing 
through the origin that does not meet this convex hull. By the same 
argument, this geometric condition is equivalent to ^ Qk-m{U)- • 

This result, combined with the motivation from quantum computing 
discussed above, naturally leads to a focus on the unitary case of the 
conjecture. We introduce the following nomenclature to delineate the 
generic unitary subcases. 
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Definition 3.4. We will use the notation Conj (A^, k) to denote the 
sub- conjecture of Conjecture A given by the statement Eq. UJ}^) for 
a given pair [N, k] and all unitary operators on N -dimensional Hilbert 
space with non- degenerate spectrum. Further, we will say Conj (A^, k) is 
constructively verified for a given pair [N, k] if it is shown that Eq. / flOj) 
holds for every unitary operator U onTi = , and if, whenever Ak{U) 
is nonempty, for every A G Ak{U) a rank-k projection P can be explic- 
itly constructed such that PUP = XP. 
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Figure 1. Conj (A^, A;) for nondegenerate unitary U. 
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4. The Structure of 

In this section we analyse the geometric structure of the set flk{U) 
for a generic unitary U on iV-dimensional Hilbert space with non- 
degenerate spectrum. Let us first estabhsh notation we will use for 
the rest of the paper. 

We shall consider the case of a unitary U with eigenvalues Xj = 
exp{i9j), j = 1,...,N, such that < < ^2 < • • • < ^'at < 27r. 
Thus, the eigenvalues Xj are ordered counterclockwise around the unit 
circle 9© in C (we use D to denote the closed unit disc). For multiple 
eigenvalues the numbering is arbitrary, but we choose an orthonormal 
system of eigenvectors {ipj) G H such that 

(11) U\^P,) = A,#,). 

When appropriate we extend the numbering of the Xj and lipj) cycli- 
cally: for example, Aat+i means Ai. Given integers i,j with i < j < 
i + N, let D{i,j, U) denote the compact convex subset of C bounded 
by the line segment from Aj to A^ and the counterclockwise circular arc 
from Xj to Aj; recall our conventions about cyclical numbering of the 
Xj. We interpret D{i,i + A^, U) as {Aj}. 

We first show that the form of VLi^{U) is simpler than what Defini- 
tion [XD suggests. In particular, Vt^iU) is a filled convex polygon with 
at most sides. 

Lemma 4.1. For all k > 1 and every unitary U G C{C^), the set 
0,k{U) is the convex polygon given by 

N 

(12) Qk{U) = f]D{i,t + k,U). 

i=l 

Proof. Let Q denote the intersection of Eq. (fT2l) . The cardinality of 
the set 

S = {l,2,...,N}\{i + l,i + 2,...,i + k-l} 

(where integers are interpreted modulo A^) is IS"! = A^ — A; + l. Thus 

ilkiU) C conv({Aj- : j G S}) C D{i, i + k,U). 

Since this holds for all i, we have flk{U) C Q. 

On the other hand, if 15"! = — A; + 1 = s we may write S = 
{ii,i2, ■ ■ ■ ,is} with 1 < ii < i2 < ■ ■ ■ < is < N and 

s 

conv({Aj : j G S}) = f^D{ij,ij+i,U) 
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(it is understood here that is+i = ii)- Since between ij and ij+i there 
are ij+i—ij — l integers that are omitted from S, we must have ij+i—ij < 
k, so that D{ij, ij+i, U) 3 D{ij,ij + k, U). It follows that fl is contained 
in conv({Aj : j G S}) for each such S. Hence flk{U) ^ fl, and equality 
is verified. ■ 

The following containments are direct consequences of this result. 
Corollary 4.2. For all U and all k, we have Qk+i{U) C Q]^[U). 

Proof. In view of Lemma 14.11 we need only observe that, for any i, 
D{i,i + k + l,U) C D{i,i + k,U). ■ 

Corollary 4.3. IfV is unitary and spec(U) C spec(y), then ^k{U) ^ 
Qk{V) for allk. 

Proof. It is enough to treat the case where spec(V) = spec([/) U 
{^AT+i} and to arrange the geometric ordering so that spec(f/) = 
{Ai, A2, . . . , Atv} and spec(l^) = {Ai, A2, . . . , Aat+i}. For 1 < i < N — k 
we have D{i, i + k,V) = D{i, i + k,U)] for N — k<i<N,we have 
D{i, i + k, V) 3 D{i, i + k, U) (note that the extended numbering within 
spec(\/) is done modulo + 1); finally, D{N + 1,N + 1 + k,V) 3 
D{N, N + k, U). Using Lemmailll we see that VtkiV) 3 Vtk{U). ■ 

The following corollary points out that Conj (A^, k) is easy to verify 
when k divides N . 

Corollary 4.4. Suppose that k divides N , with m = N/k. Then Qk{U) 
is the intersection of k m-gons, and Conj(N, k ) follows. 

Proof. Let Si = {i, i + k, i + 2k, . . . , i + {m — l)k} {i = 1, 2, ... , k); these 
index sets partition {1, 2, . . . , N}. In view of Lemma 4.1, the m-gons 

G,= f] D{i,i + k,U) 

intersect to form VtkiJJ)- Now consider A G VLkiU). Since A G for 
each i we may write 

A = ^ tijXj 

as a convex combination {tij > 0, Ylij^Si ~ -^^^ ^ ~ 2, . . . , A;, 
let 
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Figure 2. Corollary 4.4 in action: (a) ^3(U) as inter- 
section of 3 quadrilaterals when = 12; (b) Q^lU) as 
intersection of 4 triangles when = 12. 

clearly the are unit vectors, and they are orthogonal because the 
Si are disjoint. We see that (f/ — A/)|0j) ± for all i, i', so that V = 
span{|0j) : i = 1,2, . . . , k} satisfies (f/ — A/)V -L V. We can thus define 
the rank-fc projection P = Yli=i \4>i){4>i\ onto V. Then P{U — XI)P = 
and it follows that A G Ak{U). Together with Proposition 13.21 we have 
verified Conj ( A^, k) in these cases. ■ 
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Remark 4.5. The corollary above is usually of interest when m > 3. 
If N = k (m = 1), then Qk{U) will be empty unless U is a scalar. If 
N = 2k (m = 2), then D{i, i + k,U) n D{i + k,i + N, U) is the line 
segment [Aj, Aj+fc]. Thus Qk{U) is usually empty in this case for k > 2; 
in cases where spec(f/) has a special symmetry, QkiU) may be a single 
point. 

Under certain conditions Corollary 4.3 may be made more precise. 

Proposition 4.6. Suppose V is unitary on N + 1- dimensional Hilbert 
space and spec(y) = {Ai, A2, . . . , Aat+i} with distinct Xj. IfQ^+iiV) 7^ 
0, then 

N+l 

(13) Qk{V) = U QkiVj), 

where Vj is the unitary with spectrum spec(V") \ {Aj}. 

Proof. Arguing as in the proof of Corollary 4.3 we see that ^k{VN+i) = 
nk{V)nC{N + 1), where 

N 

C{N + 1)= fl Dij,j + k + l,V). 

Let A be a point in Vlk+iiV) (which is nonempty by hypothesis). Then 
A G C{N + 1) by Lemma Sil Note also that C{N + 1) includes the 
counterclockwise arc A{N + 1) of from AAr+fc+i = Xk to AAr+i_fc; 
since N + 1 — k > k (otherwise we would have N + 1 < 2k < 2{k + 1) 
so that Vtk^iiy) would be empty), the arc A{N + 1) has nonempty 
interior (relative to (9D). Likewise, for each j = 1, 2, . . . , A^-|- 1 we have 
^fc(^) — ^k{y) n C{i) where C(j) is a convex set containing A and 
an arc A{j). It remains to show that \JjC{j) = ©. The arcs A{j) 
(overlapping in general) cover all of and each C{j) includes the 
"sector" conv({A} U A{j)); these sectors certainly cover D. ■ 

The previous result raises the general question of which Qk{U) are 
nonempty. This question is of course important for the identification 
and construction of error-correcting codes. In particular, Conj ( A^, k) 
may include cases where both Ak{U) and flk{U) are empty, but this is 
not helpful for the applications. The question is somewhat clarified by 
the following. 

Theorem 4.7. Let U G £(C^) be a unitary with N distinct eigenvalues 
and let k > 1. Then we have the following conditions on Q^iU). 
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(1) IfN< 2k, then n^iU) = 0. 

(2) If N = 2k, then Qk{U) is empty if the line segments [Xj,Xj+k] 
do not intersect, and, otherwise, it is the singleton set given by the 
intersection point of these line segments. 

(3) If 2k < N < 3k — 2, then Qk{U) can be either empty or non- 
empty. 

For any unitary U G C{C^), whether or not the eigenvalues are 
distinct, we have: 

(4) If N > 3k — 2, then Qk{U) is always nonempty. 

Proof. To see (1) note that D{j,j + k, U)nD{j + k,j + 2k, U) is the sin- 
gleton {Xj+k} and that the {Xj+k} are distinct. See also Proposition 1 
from [5j. The case (2) is a special case of Corollary 14.41 Here, flk{U) 
is non-empty (and a singleton set) precisely when the line segments 
[Xj, Xj+k] intersect in a common point. 

In view of Corollary 14. 3[ it is sufficient for the statement (3) to 
provide an example of U with = 3{k — 1) and ^lk{U) = 0. Such an 
example is provided by grouping A; — 1 of the Xj close to and on the 
counterclockwise side of each of the cube roots of unity. It is then clear 
that 

D{k~l,2k-l,U)r] D{2k -2,3k- 2, U) n D{3k -3,Ak- 3, U) = 0; 

note that here = 3A; — 3 so that 3/c — 2 = 1 and Ak — 3 = k. Again 
recalling Lemma W7[\ we see that a fortiori flk{U) = 0. 

Concerning the statement (4), we recall Helly's Theorem from convex 
analysis (see, for example. Chapter 3 of |13]): a family of compact, 
convex subsets of has nonempty intersection provided any subfamily 
of size d + 1 has nonempty intersection. Since each D{i,i + k,U) is 
compact and convex in = C we need only prove that A^ > 3A; — 2 
implies 

D{a, a + k,U)n D{b, b + k,U)n D{c, c + k,U)^%, 

then invoke Helly's Theorem for d = 2 and Lemma 14. 1[ Such a triple 
intersection could only be empty if the complements in D covered all of 
D. In particular, the omitted arcs (strictly between and Xa+ki etc) 
would cover 9D. Each of these arcs contains k — 1 points from spec([/), 
so we would have the contradiction N <3{k — 1). ■ 
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Figure 3. Illustrating Theorem 4.7: (a) Here N ^ 9 
and k — 4; since N < 3k — 2 it may happen, as here, 
that n^{U) = 0; note that D{3,7,U), 0(6, 1,U), and 
D{9,4,U) overlap pairwise but have no common point, 
(b) Here N = 7 and k = 3; since N > 3k — 2, the core 
^3{U) of this 7-pointed star must be nonempty. 
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5. Verification of the Conjecture for N >3k and 
Construction of Codes for Binary Unitary Channels 

In this section we give a constructive verification of Conj(A^, k) in the 
cases N > 3k. In such cases we explicitly construct error-correcting 
codes. We require more notation. 

Let Afc(C/) denote the set of those A e C such that for some k disjoint 
subsets Si, S2, ■ ■ ■ , Sk of {1,2,..., A^} wc have A G conv({Aj : j G Si}) 
for each i. The A in the notation Ak{U) is to recall the "disjoint" 
condition. Evidently Ak{U) C Ak{V) whenever spec(f/) C spec(V^). 

A straightforward generalization of the construction used in Corol- 
lary 4.4 yields the following result. 

Lemma 5.1. For every unitary U G C{C^) and k > 1 we have 
Ak{U)CAk{U). 

Proof. Let A G Ak{U) be expressed as a convex combination of the 
eigenvalues for U for each i: 

A = ^ tijXj, tij > 0, ^ tij = 1. 

Let = v^lV'j)- Then |0i), . . . , are orthonormal (or- 

thogonal because the Si are disjoint). Let P be the orthogonal rank-A; 
projection onto V = span{|(^i), . . . , Then for each \'4>) G V, 

the vector [U — XI)\il)) is orthogonal to V. Indeed, \i i ^ i' then 
(U — A/)|0i) G span{|-^j) : j G S.i\, which in turn is orthogonal to 
because G span{|'?/'j) : j G 5*^/} and fl = 0. On the other 
hand, {{U — XI)(f)i\(f)i) = {U4>i\4>i) — A and 

{U(pi\(pi) = (^Xj^/Ujipjl y/Ujipj) = ^tijXj = A. 
jeSi j£Si jeSi 

Hence PUP = AP; that is, A belongs to Ak{U). ■ 

The converse inclusion holds in a wide variety of cases. 

Theorem 5.2. If N > 3k, then QkiU) = Ak{U). Hence, Conj(N,k) 
holds whenever N >3k. 

Proof. Let T{a,b,c) denote the "eigentriangle" conv({Aa, A;,, Ac}). We 
shall say T(a, b, c) and T(a', b', c') are "disjoint" if the index sets {a, b, c} 
and {a', b', c'} are disjoint; of course, the eigentriangles themselves may 
very well overlap. We show that N > 3k implies that every A G fik{U) 
lies in k (pairwise) disjoint eigentriangles, and so A G Ak{U) C Ak{U) 
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by Lemma [5.11 This is clear for = 1 since any convex polygon (here 
Qi{U), that is conv({Aj : j = 1,2,..., N}) is (in many ways) a union 
of triangles formed from the vertices. 

For A; > 1 we proceed by induction. The wedge W = D{1, k + l,U)n 
D{N — A; + 1, 1, U)) contains Qk{U)- Consider the eigentriangles 

T(1,A; + 1,2A;+ 1), ... , T(l, A; + (A^ - 3A; + 1), 2A; + (A^ - 3A; + 1)); 

note that (because > 3A;) 2A;-|-1 < N — k + 1 and that the last 
eigentriangle in this list is T{1,N — 2k + 1, N — k + 1). Thus the 
union of these overlapping eigentriangles covers the part of W that 
contains Qk{U), and hence covers fifc(f/) itself. Given any A G Qk{U), 
choose one of the eigentriangles from this list that contains A. Let the 
chosen eigentriangle be T(l, 6, c); note that b>k + l, c = b + k, and 
c < — A; + 1. We claim that A is also in Qk_i{W), where spec{W) = 
{Xj : j 7^ l,b,c}. Assuming this claim is correct for the moment, we see 
that the inductive step is achieved, since | spec(iy)| = A^ — 3 > 3(A; — 1) 
so that A lies in A; — 1 disjoint eigentriangles drawn from spec(iy) as 
well as in T(l, b, c), and hence that A G Ak{U). 

To verify the claim keep Lemma 14.11 in mind and note that the sets 
D{j,j + A; — 1, W) that intersect to form Qk-i{W) strictly contain one 
of the D{i,i + k,U) unless the arc omitted from D{j,j + k — 1,W) 
includes one of Ai,A6, Ac, in which case D{j,j + k — l,W) coincides 
with one of the D{i,i + k,U). The key point is that the arc cannot 
contain more than one of Ai, A^, Ac since these are separated by at least 
A; — 1 points in spec(W^). Thus, in fact, 3 Qk{U), and this 

completes the proof. ■ 

Let us discuss the construction of codes. The cases in which k di- 
vides A^ are perhaps the simplest cases in which codes can be explicitly 
constructed, as described in the proof of Corollary 14.41 The proof of 
Conj ( A^, A;) for N > 3k given in Theorem 15.21 is also constructive in the 
sense that the A G Ak{U) and the corresponding projections P may be 
found explicitly by an algorithm based on the proof. We state this in 
terms of the general binary unitary channel error correction problem. 

Let U he a. unitary on and let A; be a positive integer such 
that A^ > 3A;. Then by Theorem 14.71 (4) and Theorem 15. 2[ we have 
^kiU) = A/c(f/) = 0,k{U) and this set is nonempty. Following the proof 
of Theorem 15.21 (and recalling our earlier notation), a A:-dimensional 
correctable code for any channel of the form S = {y^/, a/1 — pU} can 
be constructed by: 
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(z) Compute flk{U) from the eigenvalues {Ai,...,AAr} of U. This 
can be done by using Lemma 14.11 in general, or by simpler means in 
special cases, such as that of Corollary 14.41 

(ii) Choose A G Qk{U). By Theorem l5.2l we can find fc-eigentriangles 
T{aj, bj, Cj), j = 1, . . . ,k, such that each contains A and there are no 
repeats in the set {aj, bj, Cj}^^^ C {1, . . . , A^}. 

(iii) As A belongs to each of the convex hulls convjAa^, A;,^., Ac^.}, for 
j = l,...,k, we can compute tij,t2j,t^j > 0, tjj = 1 such that 

A = tijXa- + t2jXbj + tsjXcy 

(iv) For j = 1, . . . , fc, put Sij = y/t^ and define (orthonormal) states 



\(f)j) = SijliJa^) + S2j\A,) + S3j#c,)- Let P = Ej=il0jX0jl- Then 



PC^ = span{|0i), . . . , \4>n)} is a fc-dimensional correctable code for £. 

Remark 5.3. It is instructive to rephrase this construction in terms 
of the number of qubits, assuming that the entire system has dimension 
N = 2"'. The codes constructed above work for k = [A^/3], which for 
n > 2 is not smaller than N/A = 2""^. Therefore these error- correcting 
codes support n — 2 logical qubits. Our construction shows that such 
codes are parametrized by the complex numbers A G Kk{U), which is 
a nonempty set explicitly determined by the eigenvalues for U as we 
have shown. This result is optimal for n > 3 in the sense that for a 
generic unitary one cannot obtain a code preserving n — 1 qubits. This 
observation follows from the fact that An/2{U) is not empty only in a 
very specific situation; if the N/2 lines joining opposite (with respect 
to the ordering number of the phase) eigenvalues ofU cross in a single 
point. 

This result considered in the context of qudits (which denote d-level 
systems) has the following implication: If N = and d > 3 then 
k = [N/3] > N/d, so the constructed code supports d — 1 qudits. 

6. Non-Constructive Verifications of the Conjecture 

In this section we derive a non-constructive verification of Conj ( A^, k) 
in the case N = 5, k = 2, and then we extend the proof to a variety of 
cases. For convenience we consider only U with distinct eigenvalues. 

To move beyond the limitations of Ak{U) we introduce Sfc(f/) as the 
set of all A G conv(spec([/)) such that for some single convex combina- 
tion 
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we have aij G C {i = 1,2, ... ,k; j = 1,2, ... , N) with the following 
properties: \aij\ — ■yjt'j for all i and j, and whenever i ^ i' we have 
both Ylij OiijC^ = and XjaijOi/j — 0. The symbol S in TikiU) is 
chosen to recall that we use a single convex combination to represent 
each A. 

It is easy to see that C Ak{U): consider any A e ^k{U) and 

(using again the above notation) let = aij\ipj). The conditions 
on the aij directly imply that the are orthonormal, that {U(f)i\(f>i) — 
A for each i, and that {U(l)i\d>i') = whenever i ^ %' . Let P be the rank- 
k orthogonal projection onto the subspace V = span{|0i), . . . , 
Clearly {U - A/)V is orthogonal to V, so that PUP = XP and A e 

Let a denote spec(t/) considered as a vector (Ai, A2, . . . , Xn) G 
In terms of a the following lemma provides a recipe for making elements 
of E2(C/). 

Lemma 6.1. Given p e C"^ such that 7^ p ± {l,a,a}, set s{p) — 
\Pj\ arid 

Then f(p)eMU). 

Proof. Let aij = \/\Pj\/s{p) and a2j = aijPj/\pj\ (if Pj = 0, let a2j = 
aij (= 0)). With A = /(p) and tj = \pj\/s{p) we easily verify the re- 
quirements for the aij. For example, Xjaija2j — (Ej \Pj)/^{p) — 

{p,w)/s{p)^o. m 

Corollary 6.2. For any k >2 we have Ak{U) C S2(t/). 

Proof. Consider A G Ak{U). Using again the notation from the proof 
of Lemma 5.1, let tu = (^11,^12, • • • ,^ijv) with the understanding that 
tij = if j ^ ^i; similarly define t2*- Let p = ti* — t2*- Then p 7^ 
since the supports of ti* and t2* are disjoint; (p, 1) = Ej ^ij — J2j — 
1 - 1 = 0; (p, cr) = ^jtijXj - Y^jt2jXj = A - A = 0; similarly p is 
perpendicular to a. Thus /(p) G T,2{U). Finally, since the supports of 
ti^ and ^2* are disjoint, /(p) = {Y.jtijXj+Y^^t2jXj) /2 = (A+A)/2 = A. 

■ 

Remcirk 6.3. Although Corollary 6.2 is sufficient for our present needs, 
it should be noted that in fact Afe(C/) C Efc(C/) for any k. Indeed, for 
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A G ^k{U) there are tij > (i = 1,2,..., A;; j = 1,2,...,N) such 
that Ylj^ij ~ ^> Ylj^ij-^j ~ ^> ^'^^ /^'^ each j at most one of the tij 
is nonzero. Thus, setting s^j = ^Jtij we have a k x N matrix S such 
that 5* diag(Aj)S'* = Xlk and SS* = Ik- Let F be any k x k unitary 
matrix with \fij\ = l/\/k (for all i,j); for example, F could be the "fi- 
nite Fourier transform", where fij = oo'^^ /y/k with cu a primitive k-th 
root of unity. Setting R = FS we have Rdia.g{Xj)R* = FXI^F* = XI^, 
and RR* = Ik- Thus we may verify that X G T,k{U) by considering 
(for each j) \aij\ is then independent of i; in fact, for 
the given j just one Smj ^ so that \aij\ = \fimSmj\ = \smj\/Vk- 

It is clear that when = 5 only the boundary dfl2{U) is captured by 
A2{U); each point on an edge of the pentagon fl2{U) belongs both to a 
line segment conv({Aj, Aj+2}) and to the complementary eigentriangle 
T{i — l,z + l,i + 3). To capture the interior of ^2{U) we turn to T,2{U)- 
The topological techniques used in the proof of the following result 
are interesting mathematically but make the construction of specific 
projections P corresponding to A G A.2{U) more difficult than we saw 
with techniques based on Ak{U). 

Theorem 6.4. If N = 5, then n2{U) = ^2{U). Thus, Conj(5,2) is 
correct. 

Proof. By Corollary [O we know that dn2{U) = A2{U) C S2(f/); to 
capture the interior we first elaborate the ideas in the proof of that 
corollary. Let Oj be the vertex of Q2{U) at the point of intersection be- 
tween the line segment [Aj, Ai+2] (in other words, conv({A.j, Aj+2})) and 
the line segment [Aj+i,Aj+3]. Each a G [ai,ai+i] also lies in [Aj+i,Ai+3] 
and so has a unique representation as a convex combination tij{a)Xj 
with tij{a) = when j ^ {i + l, i + 3}. Likewise a G [at, a^+i] also lies in 
the eigentriangle T{i, i + 2, i+i) and so has a unique representation as a 
convex combination Sij{a)Xj with Sjj(a) = when j G {i + l,i + 3}. 

For each a G [oj, aj+i], let Pi{a) = tj=K(a) — Sj*(a); as in the proof of the 
last corollary, we see that Pi{a) G X where X = {l,a, a}-'-\{0} and that 
a = f{pi{a)). Note that X = \ {0} = \ {0}, so that X is simply 
connected. Moreover, each pi is continuous on [a^, Oj+i]. Because of the 
uniqueness of the representations as convex combinations, tj*(aj+i) = 
S(j+i)*(ai+i) and Si*(aj+i) = t(i+i)*(ai+i). Thus Pi{ai+i) = -pi+i(aj+i). 
Let 7o : [0,10] dQ2{U) be a (continuous) path traversing dQ2{U) 
twice in the counterclockwise direction, beginning and ending at ai 
and such that 7o([j, j + 1]) = [0^+1,0^+2] {j = 0,1,..., 9) with the 
understanding that the are numbered modulo 5 (= A^). For t G 



HIGHER-RANK NUMERICAL RANGES 



19 



+ 1] let To(t) = {—iypj^i{jo(t)). The alternating signs ensure 
that To is continuous (even at the integers), and the double circuit 
ensures that ro(0) = ro(lO); in other words, that Tq is a loop in X. 
Furthermore f{To{t)) = /(±p,+i(7o(t))) = ToW- 

Thus, given any A in the interior of Q2{U), the winding number of 
/ o Fq relative to A is 2; that is, wx{f o Fq) = 2. Since X is simply 
connected, the loop Fq is part of a continuous family of loops F^ (0 < 
s < 1) in X such that Fi is the constant loop at some G X. Suppose 
A does not lie on any of the loops / o F^; then w\{f oTs) = 2 for all 
s, a contradiction, since w\{f o Fi) = (/ o Fi is the constant loop at 
f{p*)). Thus, for some s,t, f(Ts(t)) = A. Since Ts(t) G X, Lemma \6A] 
implies that A G T,2{U). ■ 

This result has a number of consequences, as follows. 
Corollary 6.5. For every natural number m, Conj(5m,2m) is valid. 

Proof. Let Sj = {j,j + m,j + 2m, j + 3m, j + 4m} (j = 1,2, ... ,m). 
The Sj partition {1,2,..., 5m} into m disjoint subsets. In view of 
Lemma [4. 11 we have 

5m 

^^2m(f/) = []D{z,i + 2m,U) = 

i=l 

m m 5 

f|(f| D(2,2 + 2m,[/)) = f]{f]D{z,z + 2,U,)), 

j=l i€Sj j=l i=l 

where Uj is the unitary with spectrum spec(f/j) = {Aj : i G Sj}. Thus 
^2m{U) = rijli ^2{Uj) and this is fljli ^2{Uj) by Theorem 16.41 (since 
I spec(f/j)| = 5). 

Consider any A G Q2m{U)', for each j = 1,2, ... ,m this A G A2(f/j) 
and we have a 2-dimensional subspace Vj of span({ : i G Sj}) such 
that {U — XI)Vj -L Vj. Now Vj and {U — \I)Vj are subspaces of the 
mutually orthogonal 

span{\i;i) : i E Sj} (j = 1, 2, . . . , m). 

Thus V = Vi + V2 + ■ ■ ■ + Vm is a 2m-dimensional subspace such that 
(^7- A/)V± V, so that A G A2m(?7). ■ 

Remark 6.6. Although Corollary 6.4 ensures that Conj(15,6) is cor- 
rect, for instance, it may not he useful in applications because it can 
happen thatVL^{U) = when N = 15 (15 = 3-6— 3; recall Theorem 4-6). 
Moreover, ^17{V) may be empty in dimension 16 (IQ < 3-7 — 3) so that 
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we cannot use Proposition 4-5 to support the "induction" Conj(15,6) 
=^ Conj(16,6). Thus Conj(16,6) remains undecided at the moment. 

Corollary 6.7. Conj(3k — l,k) holds for all k > 1. 

Proof. The case A; = 1 is trivial (for all A^). The case k = 2, N = 5 was 
proved in Theorem 16. 4[ Based on this we can make an induction on k 
somewhat similar to that used in Theorem 5.2. For k + 1 > 2 consider 
U with I spec([/)| = iV = 3(A; + 1) - 1. Let Wi denote the wedge 

Wi = D{1, l + {k + 1), U) n D{1 + {2k + 1), 1 + {2k + l) + {k + 1), U). 

Note that Qk+i{U) C Wi and in fact flk+i{U) is contained in the 
eigentriangle Ti = T(l, 1 + (A; + 1), 1 + (2fc + 1)), since (1 + {2k + 1)) - 
(1 + (A; + 1)) = A; < A; + 1 and 1 + (2A; + 1) + (A; + 1) = + 1 = 1. Let 
Ui he the unitary with spectrum 

spec(f/i) = spec{U) \ {Ai, Ai+(fc+i), Ai+(2fc+i)}. 

Then fifc(f/i) D l^fc+i(t/)) n D{1), where D{1) = D{l + k,l + {2k + 
2),U). Thus any A G flk+i{U) fl D{1) is in Qk{Ui) as well as in Ti. 
Since | spec(f/i)| = 3k — 1 the inductive hypothesis ensures that A G 
Ak{Ui) and there is a A;-dimensional subspace Vi of spa.'n{{\ipi) : A,; G 
spec([/i)}) such that {U — A/)Vi -L Vi. Since A G Ti we also have 
a 1-dimensional subspace V[ of spa.n{\ipi) , \ipi+{2k+i))} such 

that {U — A/)V( -L V'l- The (orthogonal) sum Vi + V[ shows that 
A G Ak+i{U). 

Similarly we have D{k+3) = D{k+3 + k, k+3 + {2k+2), U) such that 
A G ^]fc+l(f/)^D(A;+3) implies A G Afc+i(f/). Finally D(l)UL>(A;+3) = 
D, since 1 + (2A; + 2) = A; + 3 + A; and A; + 3 + (2A; + 2) < iV + 1 + A; (if 
and only if A; + 1 > 2). ■ 



Corollary 6.8. For N = 7, we at least have dn3{U) C A3{U). 

Proof. Along the lines of the proofs above, we need only show that 
A G dn3{U) implies A G ni{U') n n2{U"), where the spectra a', a" of 
U',U" partition the spectrum cr of U, \cr'\ = 2, and \a"\ = 5. Then 
A G Ai{U') and, using Theorem 6.3 we also have A G A2{U"). To 
complete the argument note that, if A belongs to one of the line seg- 
ments forming a side of dfl3{U), then A G [Aj, Aj+s] for some i and 
we set a' = {Ai,Aj+3}. Since Aj+i and Ai+2 are the only points of 
a" = a \ cr' in the counterclockwise arc from Aj to Aj+a, we must also 
have A G Q2{U"). ■ 
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Remark 6.9. // there were some sort of a priori convexity result for 
the sets Afc(T) - along the lines of the Hausdorff-Toeplitz Theorem for 
the classical numerical range (= Ai(T) ) - many of our arguments could 
be simplified and extended. For example, from Corollary 6.6 we could 
derive Conj(7,3). If Conj(N, k) is true in general then we also have 
convexity of all Ak(T) for all normal T , hut we do not know of any 
independent argument for such convexity. In fact, Conj(7,3) seems to 
be the "smallest" case that presently remains unsettled. 

7. Outlook 

As discussed above, the full Conjecture A remains open. Here we 
have constructively verified the conjecture in a wide variety of cases, 
and, curiously, non-constructively in other cases. An overarching con- 
ceptual proof covering all cases would be of great interest. A possible 
avenue to such a result could come through a general convexity theorem 
for higher-rank numerical ranges, independent of normality or unitar- 
ity, though the present work suggests that establishing such a result 
would be a delicate matter. Progress in this direction is contained in 
the recent work [4]. 

To apply the compression approach [5] to broader classes of noise 
maps (in particular to randomized unitary channels with more than two 
unitary errors) as a means to construct ideal correctable codes, a better 
understanding is required of joint solutions to the family of equations 
given by Eqs. ([2]). Furthermore, we have focussed on the generic case 
of non-degenerate spectrum to streamline the presentation. But many 
naturally arising physical examples include degenerate spectra. There 
are extra technical issues to overcome in such cases, but we expect our 
results can be extended to the case of degenerate spectra. 

It would be interesting to consider possible infinite-dimensional ex- 
tensions of these results. We have also not studied here possible impli- 
cations of the compression approach to more subtle subsystem codes 
[m, IT2j . Nor have we considered possible applications to approximate 
error-correction ^ [71 114L 116] , and specifically to noise maps that have 
two unitary errors which occur with high probability. 

Acknowledgements. M.D.C, J.A.H. and D.W.K. were partially sup- 
ported by NSERC. D.W.K. also acknowledges support from ERA, CFI, 
and OIT. K. Z. acknowledges a partial support by the grant number 
PBZ-MIN-008/P03/2003 of Pohsh Ministry of Science and Information 
Technology. 



22 M.-D. CHOI, J.A. HOLBROOK, D.W. KRIBS, K. ZYCZKOWSKI 



References 

[1] R. Alicki, and K. Lendi, Quantum dynamical semigroups and applications, 

Springer-Verlag, Berlin, (1987). 
[2] H. Barnum, E. Knill, Reversing quantum dynamics with near-optimal quantum 

and classical fidelity, J. Math. Phys. 43, 2097 (2002). 
[3] I. Bengtsson, and K. Zyczkowski, Geometry of quantum states, Cambridge 

University Press (2006). 
[4] M. D. Choi, M. Giesinger, J. A. Holbrook, and D. W. Kribs, Geometry of 

higher-rank numerical ranges, preprint, 2007. 
[5] M. D. Choi, D. W. Kribs, and K. Zyczkowski, Quantum error correcting codes 

from the compression formalism. Rep. Math. Phys. 58, 77-86 (2006). 
[6] M. D. Choi, D. W. Kribs, and K. Zyczkowski, Higher-rank numerical ranges 

and compression problems, Lin. Alg. AppL, 418, 828-839 (2006). 
[7] C. Crepeau, D. Gottesman and A. Smith, Approximate quantum error- 
correcting codes and secret sharing schemes, quant-ph/0503139. 
[8] D. R. Farenick, Matricial extensions of the numerical range: A brief survey. 

Linear and Multilinear Algebra 34, 197-211 (1993). 
[9] P. Halmos, A Hilbert space problem book, D. Van Nostrand Company, Ltd., 

Toronto, (1967). 

[10] E. Knill and R. Laflamme, A theory of quantum error- correcting codes, Phys. 

Rev. A 55, 900 (1997). 
[11] D. Kribs, R. Laflamme and D. Poulin, Unified and generalized approach to 

quantum error correction, Phys. Rev. Lett. 94, 180501 (2005). 
[12] D. W. Kribs, R. Laflamme, D. Poulin and M. Lesosky, Operator quantum error 

correction. Quantum Inf. & Comp. 6 382 (2006). 
[13] S. R. Lay, Gonvex sets and their applications, John Wiley & Sons, (1982). 
[14] D. W. Leung, M. A. Nielsen, I. L. Chuang and Y. Yamamoto, Approximate 

quantum error correction can lead to better codes, Phys. Rev. A 56, 2567 

(1997). 

[15] C.-K. Li, and N.-K. Tsing, On the kth matrix numerical range. Linear and 

Multilinear Algebra 28, 229-239 (1991). 
[16] B. Schumacher, M. D. Westmoreland, Approximate quantum error correction, 

quant-ph/0112106. 



^Department of Mathematics, University of Toronto, Toronto, On- 
tario, Canada M5S 2E4 

^Department of Mathematics and Statistics, University of Guelph, 
GuELPH, Ontario, Canada NIG 2W1 

^Institute for Quantum Computing, University of Waterloo, Wa- 
terloo, ON, CANADA N2L 3G1 

''Institute of Physics, Jagiellonian University, ul. Reymonta 4, 30- 
059 Cracow, Poland 



HIGHER-RANK NUMERICAL RANGES 



23 



^Center for Theoretical Physics, Polish Academy of Sciences, Al. Lot- 
NiKow 32/44, 02-668 Warsaw, Poland 



